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We prove the Vollhardt-Wolfle hypothesis that the irreducible vertex [/^k' il) appearing 
in the Bethe-Salpeter equation contains a diffusion pole (with the observable diffusion 
coefficient D{uj,q)) in the limit k + k' ^ 0. In the quantum kinetic equation, the 
quantity f/kk'(q) plays the role of a transition probability IVkk' and its anomalous 
growth as D{uj,q) — > is the physical reason for localization. As u; — > 0, the relation 
D{u, q) = {—iuj)d{q) holds in the localized phase, where d{q) is a regular function of g^, 
related with the properties of a typical wave function. The presence of a diffusion pole 
in ?7kk'(Q) makes it possible to represent the quantum "collision operator" L as a sum 
of a singular operator Lsmg, which has an infinite number of zero modes, and a regular 
operator Lreg of a general form. Investigation of the response of the system to a change 
in Lreg leads to a self-consistency equation, which replaces the rough Vollhardt-Wolfle 
equation. Its solution shows that D{0,q) vanishes at the transition point simultaneously 
for all q. The spatial dispersion of D{uj,q) at ^ is found to be ~ 1 in relative 
units. It is determined by the atomic scale, and it has no manifestations on the scale 
q ~ associated with the correlation length ^. The values obtained for the critical 
exponent s of the conductivity and the critical exponent u of the localization length in 
a (i-dimensional space, s = 1 (ci > 2) and z/ = (rf — 2)~^ (2 < c? < 4), u = 1/2 {d > 4), 
agree with all reliably established results. With respect to the character of the change 
in the symmetry, the Anderson transition is found to be similar to the Curie point of 
an isotropic ferromagnet with an infinite number of components. For such a magnet, 
the critical exponents are known exactly and they agree with the exponents indicated 
above. This suggests that the symmetry of the critical point has been established 
correctly and that the exponents have been determined exactly. 
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1. INTRODUCTION 



It is now widely acknowledged (see, for example, Ref. 1, p. 76) that the theory of 
phase transitions should, in principle, be constructed as a symmetry theory. Specifically, 
the effective Hamiltonian of the system is represented in the form 



where He is the critical-point Hamiltonian, possessing a high symmetry; Hint is a general 
operator which is compatible with the symmetry of the total Hamiltonian H; and, r is a 
parameter that measures the distance to the transition. The most general motivation for 
the separation (1) is that the set of Hamiltonians H^ (for example, Hamiltonians of different 
ferromagnets at the Curie point) should be separated from the set of all Hamiltonians H 
by imposing some kind of additional conditions which can be interpreted as generalized 
symmetry requirements. 

In this approach, the problem consists of determining the complete symmetry of the 
Hamiltonian H^, thus far, it has been impossible to do for most phase transitions. For exam- 
ple, the well-known Landau theory [2] starts from the obvious symmetry of the Hamiltonian 
and does not take into account scale invariance and other symmetry elements arising as a 
result of the fluctuations near the critical point (Ref. 3, Chap. 9, §2). The Landau theory 
is exact, giving an example of a complete symmetry theory, only in high-dimensional spaces 
where the additional symmetry associated with fluctuations does not arise. Another ex- 
ample is the conformal theory of phase transitions for the two-dimensional case [4], which, 
proceeding from the conformal invariance of the system at the critical point and the finite- 
ness of the number of strongly fluctuating quantities, fixes a discrete set of possible values 
for critical exponents. 

In the present paper we adopt the symmetry approach to the investigation of the An- 
derson transition [5 - 10], making a separation of the type (1) not for the Hamiltonian H 
but for an operator L which is the quantum analog of the Boltzmann collision operator. 
The theory is based on the following initial assumptions. 

1. The Schrodinger equation in a space of dimension d 



describing the motion of non-interacting electrons with an arbitrary spectrum e(p) in a 
random potential V{r) is studied. As for the random potential, it is assumed only that the 
averages with respect to its realizations can be calculated by the diagrammatic technique. 
The existence of a diagrammatic technique for the standard models of a random potential 
can be proved directly [1, 11, 12]. In the general case, the limits of applicability of the 
diagrammatic approach are poorly investigated: some problems certainly arise for quasi- 
random systems [13, 14]. 

The exact Green's function of Eq. (2) is expressed in terms of the eigenfunctions ■0s(r) 
and eigenvalues (s = 1, 2, N), 



H = H, + tH, 



int 





(2) 




(3) 
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The averaged Green's function (G{r\r')) is determined by a diagrammatic series (Fig. la), 
and in accordance with current ideas [8, 15] it is assumed to be analytic at the point of 
the Anderson transition: For d > 4 this was recently proved by the present author [16, 17]. 
The quantity 

0(rir2, r3r4) = (G'i+Jrir2)G':^(r3r4)) (4) 

which contains information about the kinetic properties, has a singularity at the transition 
point. This quantity is determined by a series of diagrams with four legs, constructed on 

and hues (Fig. lb), and its properties are similar to those of the two-paricle Green's 
function in the theory of interacting particles [11]. It satisfies the Bethe-Salpeter equation, 
containing an irreducible vertex U (Fig. Ic). 

2. The following symmetry elements are assumed: 

(a) Spatial uniformity in the mean. This leads to a conservation law for the external 
momenta in the diagrams. This makes it possible to express (G) in terms of the self-energy 
E 

- ek - Ek' 

and to introduce for the function (j) the three-momentum notation 0kk'(Q) (Fig- Id) and to 
write the Bethe-Salpeter equation (Fig. Ic) in the form 

0kk'(q) = G^^^,^Gt^,^{N5^_^> + l^C/^^^(q)0^^^,(q)}. (6) 

ki 

Here and below the energy variable is equal to £■ + a; for the functions G^ and E for the 
functions G"^. 

(b) Isotropy and inversional invariance in the mean. With this symmetry taken into 
account, the function 0kk'(q) depends only on scalar products constructed from k, k' and 
q, whence, specifically, 

0kk'(q) = 0-k-k'(-q). (7) 

Similarly, G^ and G^ depend on k^ and are even functions of k. 

(c) Time-reversal invariance. This property makes it possible to choose real cigcnfunc- 
tions ips{i^) and to drop the conjugation sign in Eq. (3). Then G{r,r') = G{r',r) and 
interchanging ri,r2 and r3,r4 in Eq. (4) gives in the momentum representation 

</'kk'(q) = 0-k' -k(-q) (8) 

0kk'(q) = 0(k-k'+q)/2,(k'-k+q)/2(k + k') (9) 

Comparing Eqs. (7) and (8), we obtain 

0kk'(q) = 0k'k(q)- (10) 

Solving Eq. (6) formally for the function ?7kk'(q) using Eqs. (7) and (10), it is easy 
prove similar properties for this function: 

C/kk'(q) = C/-k,-k'(-q), t^kk'(q) = t^k'k(q). (11) 
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Figure 1: a, b — Diagrams for the average Green's function (a) and the quantity (b), 
which correspond to a Gaussian random potential [1] or a Born approximation for randomly 
distributed impurities [11] (in what follows, their specific form is not used), c — Graphical 
representation of the Bethe-Salpeter equation; d — explanation of the three-momentum 
notation. 
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3. It is conventionally assumed that the Anderson transition occurs from a phase with 
exponential localization of the wave functions into a phase with a finite diffusion coefficient. 
The existence of exponential localization in the limit £^ — > — oo and finite diffusion for large 
positive values of E (for d > 2 and an unbounded spectrum e(p), < e(p) < oo) has been 
firmly established, as a result of many investigations, for Eq. (2). The proof of the existence 
of a mobility edge is based mainly on Mott's argument [7]: The existence of states with 
different degree of localization and the same energy is impossible because of instability with 
respect to an infinitesimal perturbation of general type. Mott's argument does not forbid, 
however, the existence of intermediate states — with power-law localization, hybrid states, 
and so on — and correspondingly different types of " Anderson transitions" (for example, 
in the quasi-random systems [13, 14] the transition occurs from exponential localization to 
a ballistic regime). In the present paper the first instability, arising with a motion from 
deep in an exponentially localized phase, is investigated and it is shown that it does indeed 
correspond to a transition into a phase with finite diffusion. 

4. The general notions of the modern theory of critical phenomena [1] (such as param- 
eter space, critical surface, relevant and irrelevant parameters) are used. 

5. The theory is based on the physical idea that the localization phenomenon is asso- 
ciated with a diffusion pole in the irreducible four-leg vertex 

uMa) = f/S?(q) + us'^t) = c'S?(q) + (12) 

proposed by Vollhardt and Wolfic in the so-called "self-consistent theory of localization" 
(see Ref. 18, and also Refs. 10 and 19). 

This idea agrees with the theory of weak locahzation [20 - 22], according to which 
the diffusion pole in C/kk'(q) determines the main quantum corrections to the conductivity 
which in turn determine the scaling behavior in a space with dimension d = 2 + e. The 
diffusion pole in f/kk'(Q) with the classical diffusion coefficient Dd arises as a result of 
summation of fan-shaped diagrams [20]; Vollhardt and Wolfie conjectured that when all 
diagrams are taken into account, Dd is replaced by the exact diffusion coefficient D{u,q). 
They approximated C/kk'(q) by Eq. 12 with C^kk^(q) = const, F(k, k', q) = const and solved 
approximately the Bethe-Salpeter equation (6), which, using the Ward identity [18] 

^Sk(q) = ^E^kk,(q)AGk,(q), (13) 

ki 

AGk(q) = - , AEk(q) ^ E^^^/^ - ^t^/, (14) 

was rewritten in the form 

[-a; + (ek+q/2 - ek-q/2)]0kk'(q) + 
+ ^ E t^kk'(q) [AG'ki(q)0kk'(q) - AG'k(q)0kik'(q)] = AG'k(q)A^5k-k'. (15) 
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There exists a simple estimate which gives the same results. We note that the second term 
in the left-hand side of Eq. (15) is reminiscent of a Boltzmann collision integral and indeed 
transfer into it in the limit of weak disorder (Sec. 3). It is significant that in the quantum 
region the quantity f/kk'(Q) plays the role of a "transition probability". Using analogue of 
r-approximation, D oc / oc (f/)~^ (/ is the mean free path, (...) denotes averaging over the 
momenta), and taking into account Eq. (12), we obtain the self-consistency equation of 
the VoUhardt-Wolfle theory 

D ~ const f C/o + Fo / . J^, . (16) 

\ J -icu + D{uj,q)q^ J 

This estimate is no less accurate and demonstrates more clearly the crux of the matter 
than the approximate solution given in Ref. 18 for Eq. (15). As the degree of disorder 
increases, the "transition probability" increases anomalously as a result of a decrease of the 
diffusion coefficient, making it possible for this coefficient to vanish. Neglecting the spatial 
dispersion D{u!, q), Eq. (16) makes it possible to determine the critical exponents for the 
conductivity a and the localization length ^ 

a - T^ T-" (17) 

(r is the distance to the transition); setting D — const(a;) ~ cr in the metal phase and 
D ~ (— io;)^^ in the localized phase, we obtain 



s = l, d> 2 ; V = < 



1 



d-2 



1 

2 



, 2<d<A 

(18) 

d>A 



The drawbacks of the self-consistent theory [18] can already be seen from the exposition 
given above: 

(a) The method used to solve the Bethe-Salpeter equation is rough; 

(b) The spatial dispersion D{uj,q) is ignored, while it can change substantially the 
estimate of the integral in Eq. (16); 

(c) An approximation is used for ?7kk'(^) that leads to a singularity ~ 1/a; on the 
right-hand side of Ward's identity (13) in the localized phase, and this is incompatible with 
regularity of E at the transition point. 

One of the most interesting questions in the theory of localization is connected with the 
drawback (b). It follows from the Berezinskii-Gor'kov criterion [24] that in the localized 
phase D{0,q) = (Sec. 4). The question arises, how the spatial dispersion of D changes 
near the Anderson transition. VoUhardt and Wolfle had in mind that D{0,q) vanishes at 
the transition point simultaneously for all values of q; in the clear form, such hypothesis 
was stated by Efetov [21]. The vanishing of the whole function cannot occur accidently and 
it must be supported by a deep symmetry. Does this symmetry exist? What is its nature? 
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Another fundamental question of the theory is touched upon in the drawback (c). If 
the diffusion pole in C/kk'(Q) exists, then why is there no l/o; singularity on the right-hand 
side of Ward's identity (13)? The condition for this pole to be cancelled imposes stringent 
requirements on the approximation employed, while the satisfaction of Ward's identity has 
actually never been checked in any of the existing theories [10]. 

A theory free of the drawbacks (a-c) and answering the questions raised is expounded 
below. The first half of this paper follows the scheme of the VoUhardt-Wolfle theory and 
contains a proof of the relations obtained in Ref. 18 by a chain of hypotheses or doubtful 
approximations. In Sees. 2 and 3 the diffusion poles of </>kk'(c[) ^-^e separated out and 
the result (12) is proved. In Sec. 4 the general properties of the diffusion coefficient and 
its relation with the localization of the wave functions are determined. The content of 
Sees. 5-7 replaces the rough solution of the Bethe-Salpeter equation [18]: In Sec. 5 a 
hierarchical structure of the spectrum is obtained for the quantum collision operator L; a 
separation of the type (1), convenient for a symmetry analysis, is established; a condition 
on the transition point is found and a self-consistency equation, replacing Eq. (16), is 
derived. The diffusion coefficient D{uj,q) is sought with no assumptions on the character 
of the spatial dispersion (Sees. 6 and 7), but only the solution with a weak q dependence 
is found to be internally consistent. Such weak dependence does not affect the estimate of 
the integral in Eq. (16) and leads to the result (18) for the critical exponents. In summary, 
all basic results of Ref. 18 are found to be correct, which is surprising for such a rough 
theory. 

The theory expounded starts from the obvious symmetry of the system and the ad- 
ditional symmetry of the critical point is determined in the course of the analysis. The 
inevitable question is: Are the hidden symmetry elements completely determined? There 
is a serious argument, that the determination is complete (Sec. 8): with respect to the 
character of the change in symmetry, the Anderson transition is found to be similar to the 
Curie point for an isotropic n-component ferromagnet in the limit n — > oo. This model of a 
ferromagnet is the basis of the 1/n expansion [1], its critical exponents are known exactly, 
and they are in exact agreement with Eq. (18). The isotropy of the equivalent ferromag- 
net is the symmetry that makes D{0,q) vanish simultaneously for all q. The approximate 
(to an accuracy ~ u) orthogonality of the singular part of ^/kkK^l) ^Gk(q) results in 
cancellation of singularity in the right-hand side of the Ward identity [Sec. 5.3, Eq. (13)]. 

Another method for checking the completeness of the symmetry found is to compare 
with the results of model investigations. The hypothesis that the exponents (18) are exact 
was actually stated in Ref. 25 on the basis of an analysis of all known results: 

(a) For d = 2 + e Wegner's relation s = {d — 2)i/, following from the existence of 
one-parameter scaling [21], is valid and the e expansion for the exponent i> has the form 
[26] 

i/= - + 0-e° + 0-e^ + O(e^), (19) 

which agrees with Eq. (18), if the coefficients of the higher order powers of e are also zero. 

(b) The result (18) separates the dimensions of the space dd = 2 and dc2 = 4, which 
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on the basis of independent analysis are considered to be the lower [21] and the upper (see 
the discussion and references in Ref. 16) critical dimensions. 

(c) The entire experience of the theory of phase transitions shows that for d > dc2 the 
critical exponents do not depend on d, which is the case in Eq. (18). 

(d) The exponents (18) agree with the results for d = (yo : u = 1/2 (Refs. 27 and 28) 
and s = 1 (Ref. 29); the disagreement with the result s = oo, obtained in Ref. 28, is 
discussed in Sec. 9. 

The value u — 1 oi the exponent ior d — 3 agrees satisfactorily with the results of 
numerical calculations (z/ = 1.2 ± 0.3 (Ref. 30) and u = 1.5 ± 0.2 (Ref. 31) and the 
qualitatively behavior of z/ as a function of d agrees with the estimates from hierarchical 
models [14]. In Wegner's work [32] a finite contribution ~ is obtained in Eq. (19). 
This makes the agreement with the numerical calculations of Refs. 30 and 31 much worse. 
However, this result was derived for the zero-component cr-model, whose correspondence 
with the initial disordered system is controversial (Wegner himself [32] does not reject this), 
and it is apparently correct only in the lowest orders in e (Sec. 9). 

A qualitative result of this work, which can be checked experimentally, is the assertion 
that there is no spatial dispersion of D(a;, q) on the scale in the limit a; — > (compare 
with Refs. 33, 34, and 35). The absence of significant spatial dispersion does not contradict 
the strong dependence of the diffusion coefficient Dl of a finite system on its size L [21]. 
This dependence is connected with the temporal dispersion and is determined in terms of 
the known function D{u!, q) by the relation Dl ~ D^Dl/L"^, 0) [18]. 



2. RELATIONS BETWEEN THE QUANTITIES IN THE PRESENCE OF 
SPATIAL AND TEMPORAL DISPERSION 



In this section the existence of a diffusion pole in the quantity 

1 



'^(q) = i^E'^kk'(q), (20) 

kk' 



which is the Fourier transform of the quantity (4) for coinciding arguments ri = r4, r2 = r^, 
will be proved. In contrast to Ref. 18 and other works, it will not be assumed that q is 
small. In view of the great confusion in the literature, we shall give a complete summary 
of the formulas which are relevant here. 

We shall consider the response of a system to an electric field E(r, t) ~ e*'i''"~*'^*. The 
frequency lu is assumed to be finite only in order to remove the uncertainties which appear; 
the limit a; — > is taken in the final results. Neglecting magnetic effects, the field E is 
a purely potential field. This makes it possible to confine attention to the longitudinal 
components of the susceptibilities (Ref. 36, §103). In the presence of spatial dispersion two 
definitions of the conductivity are possible: 

j(a;,q) =a(o^,q)E(u;,q) 

Uoo, q) = aicv, q)E(a;, q) ^ " + J<^'^^' ^^'> 
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which relate E with the total current j or its electric component je; the diffusion component 
of the current j^iffi^, q) ~ — iq£>(a;, q)p(a;, q) is due to the deviation of the electron density 
p from the equilibrium density. This deviation is determined by the polarizability a (</? is 
a scalar potential): 

p{uj, q) = a{uj, q)(p{uj, q). (22) 

The conductivity a appears in Kubo's formulas (see below), which determine the total 
response of the system to the field E. The conductivity a is related with the diffusion 
coefficient D by the Einstein relation 

a{u;,ci)^e^N{eF)D{u;,q) (23) 

since the change in the scalar potential (p and the chemical potential p with p{r, t) + 
eip{r,t) — const does not destroy the thermodynamic equilibrium (N^ep) is the density 
of states at the Fermi level). A relation between a, a, and a follows from the continuity 
equation: 

—iua(u, q) = [—ico + D(uj, q)q^]a{cu, q), 

ijja{uj,q) = —iq'^a{uj,q), (24) 

so that the difference between cr and 5" is important only for q 7^ 0. Using the relation (23), 
we obtain for the polarizability a and the permittivity e 



—iu + D{uJ, q)g2 



/ N T 477 47rCT(a;,q) 
e w,q = 1 - — a a;,q = 1 . 25 

It is clear from Eqs. (23 - 25) that if the diffusion coefficient D{ijj, q) is given, then all 
quantities introduced above can be determined. 

The quantities a and a are given by Kubo's formulas (Ref. 2, §126, and Ref. 37, §75) 

^(o;, q) = - / dte'^' j dre-'^^-'Oir, t)j(0, 0) - 3(0, O)j(rt)), (26) 


1 °° 

a{u;, q) = - / dte'^' j dre-^-^(j(r, t)p(0, 0) - p(0, 0)j(r, i)), (27) 

00 

a{uj, q) = / dte'''' J (ire-*'i "(p(r, t)p{0, 0) - p(0, 0)p(r, t)) (28) 


which determine, respectively, the response of the current to a vector potential, the response 
of a current to the scalar potential, and the response of the density to a scalar potential. 
The equivalence of Eqs. (26) and (27) and the relation (24) between a and a follow from 
the continuity equation for the density operator p and the longitudinal component of the 
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current operator j and the asymptotic expressions for a and a in the hmit oo oo (Ref. 
36, §78). 

We note that according to the precise meaning of Kubo's formula (see the detail dis- 
cussion in Ref. 38), the response of the system to the field D produced by external charges 
must be calculated. In this approach the Coulomb interaction between the electrons must 
be necessarily included in the Hamiltonian to avoid contradictions in the Maxwell's equa- 
tions; Kubo's formulas have a form that is somewhat different from Refs. 26-28 (Ref. 37, 
p. 413), and the correlation functions appearing in them must be calculated taking into 
account the Coulomb interaction. A different approach [38] is more convenient: The in- 
teraction between the electrons is divided into a short-range and slowly-varying long-range 
parts; the first part is included explicity in the Hamiltonian and the second part is taken 
into account as a self-consistent field, leading to screening of the field D; for this reason, 
the response to a real physical field E is studied and the correlation functions appearing in 
Eqs. (26-28) are calculated only taking into account the short-range part of the interaction. 
The latter part can be taken into account in the spirit of the Fermi-liquid theory. We shall 
neglect it completely, since in its classical formulation the Anderson transition problem 
is a problem of noninteracting electrons. We emphasize, that the word "noninteracting" 
must be understood precisely in the sense indicated above, since otherwise the concept of 
conductivity cannot be introduced in a consistent manner. 

The correlation function in Eq. (28) for noninteracting electrons in a random potential 
is calculated similarly to the correlation function for a Fermi gas (Ref. 2, §117) using, 
instead of the plane-wave representation, a representation in terms of the eigenfunctions 
V'.(r) of Eq. (2)0: 

aiu, ^) = e'l del du' ^'^'^ ' ^fj/U (e) (p.p.w)^. (29) 
J J 00 — oo' + to 

—CO — oo 

Here /o(e) is the Fermi function, and {pEPE+uj)q, is the Fourier transform of the Berezinskii- 
Gor'kov spectral density [24] 

{pEir)pEU^')) = ]^(E rsir)M^)rA^') U^')5{E - es)6iE - e., + uj)). (30) 

For small oo and zero temperature, taking the imaginary part of Eq. (29), we obtain the 
inversion of (29) 

(p,,p,,+.,)q - -— 5-^^ (31) 

(in the absence of interaction '?/'<j(r) and do not depend on eE and ep can be replaced by 
E). The standard diffusion form for {pEPE+uj)q (Refs. 10 and 24) is obtained by substituting 
the expression (25) into Eq. (31), assuming D{oo, q) is real, which in general is not the case. 



"'^For definiteness, all quantities refer to the same spin projection. For purely potential scattering the 
spin subsystems are independent and the number of spin components can be easily taken into account in 
the final results. 
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The following expression can be easily obtained for the function 0(q) {E and u are 
parameters appearing in Eq. (4)): 



oo oo , , , , 

m =fdJ dj- — ^'''<^f-^->", — -. (32) 

^ ' J J (E + uj -e + iS)(E -uj' -e-tS) ^ ' 

— oo — oo 

The polarizability Q;(c<;,q) is a generalized susceptibility (Ref. 2, §123; Ref.36, §103) and 
the oddness of ImQ;(a;, q) as a function of the frequency makes it possible to write 

ImQ;(u;,q) = ^ '—^ , Q;(a;, q) = q;(u;, q) - q;(0, q). (33) 

Substituting the expressions (31) and (33) into Eq. (32), we obtain integrals with Q;(a;',q) 
and a (—a;', q) that converge separately. Making the substitution — > —lo' in the second of 
the integrals and shifting upwards the contour integration over lo\ and taking into account 
the fact that Q;(a;, q) is analytic in the upper half-plane, we obtain 

oo _ / II \ r\ • 

0(q) duj — ..w . ^ —aE{uj,q), 34) 

J (uj" — uj — io)ujj" + lo) e^uj 

— oo 

where uo" = e — E. The second equality follows by neglecting uj" in the argument E+ui", as a 
function of which appreciable changes in a occur on the atomic scale and are not important 
in the region uu" ~ u, which makes the main contribution to the integral. Substituting 
a{u, q) in the form (25), we obtain 

where the contribution (/>reg(q) originates from the region of large values of a;" in Eq. (34) 
and is regular in the limit cj,q 0. In the localized phase, when D{uj,q) ~ {—'i^^) (Sec. 
4), to obtain the expression (35) we should add a small real frequency-independent term in 
the denominator, in order for all expressions to be meaningful; it is essential in separating 
a{0, q) from a{u!, q). In conclusion, the quantity 0(q) has a diffusion pole which contains 
the observable diffusion coefficient. 



3. SEPARATION OF DIFFUSION POLES FROM THE 
BETHE-SALPETER EQUATION 

We introduce the operator L, which is the symmetrized version of the operator on 
the left-hand side of Eq. (15), which arises as a result of the replacement 0kk'(q) ~^ 
0kk'(q)\/^^k(q) ^-iid division of Eq. (15) by ^AG'k(q): 

L(q)=Lo(q)+M(q) (36) 
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LoV'k = 4 E t/kki(q)[AGk,(q)V^k - VAGk(q)AGk,(q)V^kJ, 



ki 



M^k = (ek+q/2 - ek-q/2)V'k- 



The operator L acts in the complex space and, by virtue of Eq. (11), it is symmetrical 
with respect to the scalar product 

('/''^) = ^E'/'kV'k. (37) 
-'^ k 

Its cigcnfunctions e^'*''(q) form a complete orthonormal basis, and the eigenvalues Xsiq) 
are, generally speaking, complex. In terms of and e'^''^ the formal solution of the Bethe- 
Salpeter equation (15) has the form 



ft\ci) = V^G'k(q)eL^)(q). (38) 

At least, one eigenvalue — for definiteness Ao(q) — behaves as Ao(q) ~ for small q. 
Indeed, the operator Lq has a zero mode ■0k(q) = Y^AGk(q) and, treating the operator 
M ~ g as a perturbation, we can construct the iterative series 

eL°)(q) = const[4°Hq) + V'^'Hq) + 
#(q) = V^Gk(q), 
Ao(q) = Ai^)(q) + A(')(q) + ..., A('^)(q), ^ g", (39) 

(^(o),^W) = 0, ny^O 

in the Brillouin-Wigner form [39]. The eigenvalues As(q) are even with respect to q (see 
Appendix) and the correction Aq^'' is equal to zero, as one can easily verify directly. To 
second order in q we have 

(V>W,M^W) + (V>W,MV>W) 

^'^'^> p^V^M) 

where ■0^''^^ satisfies the equation 

-P^Mij^^'^ = LoV'(^) (41) 

and Pj_ is a projection operator onto the space orthogonal to ■^^"^ (Ref. 39). Making the 
substitution 

rf^(q) = -V^^k(q)q-lk (42) 
12 



and noting that 

1 Y: AGk(q) = ^ E 2^ImGk = -2mN{E) (43) 



we rewrite (40) in the form (vk is the velocity of electrons with momentum k) 



Ao(q) = 



27rN{E) 



1 E(q • Vk)(q • Ik) (-0 AGk(q) + ^ E(q " Vk) AG'k(q) 

k k 



(44) 



For an isotropic spectrum e(k) = /2m the expressions (44) and (41) assume, to lowest 
order in q, the form 

Ao(q) = -iL>(0, 0)?^ L'(0, 0) = a(0, Q)e-''N-\E), (45) 
-(0, 0) = ^1 E(v. • 10^AG.(0) + ^1 E ReC.« (46) 

Vk = ^ E^^kk'(0)AGk'(0)(lk - IkO- (47) 
-'^ k' 

In the limit of weak disorder, when 

AG'k(O) ^G^-G^^ 2ilm- ^ ^ -2mS{E - ek) (48) 

-£/ - Ck + ^7 

we obtain from Eqs. (46) and (47) 

^i^^^) = ^^E^^-^^SiE-e^), (49) 

Vk = ^ E 27rC/kk'(0)(lk - lk')<^(^ - ek'), (50) 
-'^ k' 

i.e., cr(0,0) is the classical conductivity, D{0,0) is the classical diffusion coefficient, and Ik 
is the vector mean free path length, determined by the standard classical equation (50) 
for scattering by impurities [40] . The results (46) and (47) extend the concept of a kinetic 
equation and a mean free path into the quantum region. The differences from the classical 
equations reduce to the following: 

(a) The 5-function expressing the energy conservation law is smeared; 

(b) the transition probability is replaced with 27r[/kk'(0); 

(c) a acquires a quantum correction [last term in Eq. (46)] of the order of the Mott 
minimum conductivity [7]. 

It is obvious from Eq. (45) that the diffusion pole is related with the zeroth term in the 
sum in Eq. (38). To compare with Eq. (35), we sum the expression (38) over k and k': 



-a; + Ao(q) n^-^ + As(q)' 
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^(q) = ^Ey^^4^Hq)- (51) 



k 

JO), 



Neglecting in Eq. (36) the operator M, we have (q) ~ YAG'k(q), whence ^ol^) ~ 
—2iTiN{E), As{q) =0, s ^ 0. Taking M into account by perturbation theory, we obtain 

Mlf~Q\ "^0 (52) 

and comparing the expressions (51) with Eq. (35) gives 

D{u;, q)g2 = iAo(q)[l + B{q)] - iu;B{q), 

<PreM) ~ q'. (53) 

The decomposition into regular and irregular parts is not unique and admits a "gauge 
transformation" 

0rep(q) = 0rep(q) - 27rN{E)C(q), 

N 2 D{uj,c[)q'^ + iu;C{c[)[-iuj + D{uj,q)q'^] 
U[uj, q)q — - 



l + C{q)[-iu; + D(u;,q)q^] 

C{q) ~ q' (54) 
up to which the identity (53) is valid. For this reason, it is convenient to set by definition 

Ao(q) = -iL'(u;,q)g' (55) 

making the assumption that the diffusion coefficient D{u!,q) determined in this manner 
is related to be observed diffusion coefficient Doi,s{^: q) by relations of the type (53) and 
(54). For any B{q) and C(q), we have £'(0,0) = Dobs{0,0), and D{0,q) and L'ofe,(0,q) 
vanish simultaneously. In practice, the difference between D{uj,q) and Dobs{^jJ,'i) is not 
important. The point is that the spatial dispersion of D{u!,q) on the scale g ~ A (A is 
a parameter of the order of the inverse interatomic distance) is of little interest; only the 
"anomalous" dispersion, determined by the scale which can arise near the Anderson 
transition, is of interest. The quantity B{q) does not contain anomalous dispersion, since 
it is determined by the function AGk(q), which is regular at the transiion point, and the 
function e[°^(q), which can be assumed to be constant (Sec. 5.4); this is also true of the 
quantity C(q), relating, according to Eq. (54), two regular functions (see, however. Sec. 
4). On the basis of what we have said above, the expression (38) assumes the form 



•^(0)/ \fiO)( \ 
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Figure 2: a — Structure of the diagrammatic series for 0kk'(q)- b — same, with the 
upper G hne inverted; the U and C/ blocks are topologically equivalent, but they correspond 
to different values of the momenta. 

It follows from the relation (9) that 0kk'(^) contains a diffusion pole in the limit k + k' — > 0, 
which can be separated from 0kk'(^)- 

, , ) ./r(q).ff (q) , -/f.,,;.(k + k-)ir..,..(k + k-) 

'^kk'iqj - -——-FJ---T-^ + 7-7r-F^777-m:j^TrvTJV2 — + '^kk'iqj- i^O 



In diagrammatic language, the pole in the limit ^ — > is related to the fact that for 
the diagrams containing two or more blocks U (Fig. 2a), the contour of integration in the 
integrals 

is confined between the poles of two Green's functions. For small U the divergence in 
the expression (58) in the limit a;, g — > is limited only by the small damping ImE in 
the denominators of the G functions and compensates the smallness associated with the 
addition of an extra block U ; all diagrams in Fig. 2a are found to be of the same order, 
and the series diverges, leading to a diffusion pole. For arbitrary U the divergence of 
the series in the limit uj. q ^ is guaranteed by the Ward identity (13). It is important 
that the diffusion pole is determined by diagrams with a large number of blocks U. Since 
G(ri,r2) = G(r2,ri) the result for 0kk'(q) "^i^^ remain unchanged, if in constructing the 
diagrams the upper G line is reversed; then the diagram contains blocks U (Fig. 2b), 
topologically equivalent to the blocks U, but taken for other values of the momenta. Now 
the poles of the two G-functions converge toward one another as k + k' 0, giving a second 
diffusion pole in the expression (57). When the upper G line in Fig. 2a is reversed, the 
diagrams containing two or more U blocks become irreducible and enter into a U block (Fig. 
2b) and, conversely, reversion of the G line in the diagram with one U block generates the 
entire sequence of diagrams in Fig. 2b with more than two U blocks. For this reason, the 
second pole term in the expression (57) is contained, with no changes, in C/kk'(*l)' differing 
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only by the contribution of the four external G lines. The result (12) with a function 
F(k, k', q) of the form 

F(k, k', q) = ^/(^Lk.+^)/2(k + k')/(^Lk+q)/2(k + ^^){Gl^^/,Gt^,^G^>^^,^Gt,^^/,)-' (59) 
is vahd for C/kk'(q)- This proves the VoUhardt-Wolfle hypothesis. 

4. BEREZINSKII-GOR'KOV CRITERION AND ITS CONSEQUENCES 

The spectral density (30) contains a singular contribution ~ 6{uS), originating from 
terms with s = s', which is finite in the localized phase and vanishes in the delocalized 
phase in the thermodynamic limit. This is the Berezinskii-Gor'kov localization criterion 
[24]. The d{u!) singularity in {pePe+uj)^ leads to, by virtue of the Eq. (32), a l/uj singularity 
in the function 0(q) (Ref. 10) 

0(q) = ^^^^A(q) + 0.e.(q), (60) 



^(q) = j dre-^'»"'^(r). 



^^'^ ^ We)^^ |V.(r)nV.(0)r<5(£; - 6,)). (61) 

A number of important consequences follow from the relation (60). 

1. Comparing Eqs.(60) and (35) shows that in the localized phase D{uj,c\) ~ a;. A 
slower dependence would destroy the 1/uj singularity in Eq. (60) and a more rapid de- 
pendence would cause the dependence on q to vanish in the singular part (35); such a 
dependence obviously exists according to Eq. (60). This result, valid in the -D(u;, q) gauge, 
in which the functions 0reg(q) in Eqs. (35) and (60) are identical, remains valid in any 
other gauge [see Eqs. (53) and (54)]. Therefore 

D{u,q)^{-iu)d{q), (62) 

where it is assumed that the limit a; — > is taken in the function d{q). Therefore it follows 
from the Berezinskii-Gor'kov criterion that D(0, q) vanishes for all q. This completes the 
proof of all of the main localization criteria [6, 10]. In view of Eq. (62), the second diffusion 
pole in Eq. (57) leads to the singularity l/cj in the sum over s in Eq. (51). To eliminate this 
singularity from 4>reg{^) the expression (54) must be transformed with C(q) ~ 1/a;, without 
destroying the proportionality ol D{uj, q) and D{uj,q} to the frequency. The function C(q) 
is determined, by virtue of Eq. (57), by the quantities AGk(q) and e^^q), which are 
regular at the transition point and do not lead to anomalous dispersion, while the associated 
renormalization of D(a;,q) is small near the transition because of the divergence of d{q) 
(see below). 
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2. A relation between the diffusion coefficient and the properties of the wave functions 
follows from Eqs. (35), (60), and (62): 

A(q) = / dve-''^'A{r). (63) 



1 + d{q)q^ 

Exponential localization of the wave functions leads to exponential decay of ^(r) at large r 
[see Eq. (61)] and the finiteness of the coefficients in the expansion over q of the right-hand 
side of Eq. (63). Because of isotropy in the mean, there are no odd powers of q and d{q) is a 
regular function of g^. It is important that this function does not contain noninteger powers 
of q, which arise naturally in the case of diffusion over fractal structures [41]. The reality 
and positiveness of d{q) follow from the reality of ^(q) and the inequalities < ^(q) < 1 
[10, 24]. 

3. Restrictions on the form of the spatial dispersion of D{uj, q) follow from the relation 
(63). In the localized phase the spatial dispersion is determined by the expansior|§ 

1 + diq)q' = f° + f^q' + f'q'' + ... + f-q'^ + 

/3o = 0, (64) 

where /?„ > 0, since the contributions, associated with the atomic scale and corre- 
sponding to Pn = 0, obviously exist. 

Different estimates show that the smoothed (over oscillations) behavior of the squared 
modulus of a typical wave function has the form 

|^(r)p = const| p ^-'-/-^ . (65) 
[ exp(-r/4^), r ~ ^ 

This behavior should actually be expected on the basis of Poincare's theorem on the analytic 
dependence of the solution of a differential equation on a parameter. If the behavior of 
the wave function at the transition point is characterized by the exponent b {0 < b < oo), 
|^c(r)P ~ then near the transition we have ip{r) ^ ^c(i') for sufficiently small r, as 
a consequence of Poincare's theorem. The theorem is valid only for a finite region, whose 
maximum size is determined by the scale ^ on which the exponential decrease of iplr) starts. 
By virtue of Eq. (61) the function A{r) has a similar behavior 

A(r) = const <^ ' ' > ^ , (66) 

I exp(-r/^), r ~ ^ 

where const is chosen from the condition A{q) = 1 at g = 0. The series expansion of A{q) 
AicO = e + + Pg" + - + e-g'" + 70 = (67) 



^In expansions of the type (64) arbitrary coefficients are assumed. Taking tlicm into account is beyond 
accuracy of tlie present analysis. 
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Figure 3: Possible configurations of the exponents and 7„ for C < (a) and C > (b). 



and the estimate of the integrals arising in Eq. (61) show that only two variants are 
possible: (a) 7„ = 2n for ^ < and (b) 7„ = max{0, 2n — (} for C > 0. Substituting the 
expressions (64) and (67) into Eq. (63) gives a relation between 7„ and 



7„ = meiXi+j+k+...=n{Pi + Pj + Pk + ■■■}, 



(68) 



leading to the two possibilities for the exponents /?„: Pi = 2, /9„ < 2n for C < (Fig. 3a) 
and Pn = maxjO, 2n — C} fo^^ C > (Fig. 3b). For these results to be valid it is important 
only that if the integral of y4(r)r" diverges in the limit — > cxd as then the integral of 
jl(^j.yn+m should diverge as since it is determined by the region r ~ ^. The specific 

approximation (66) is actually not used, but it is convenient for interpreting the results. 

To determine the localization length ^ from the known diffusion coefficient, in general, 
it is necessary to know all the exponents Pn- The result D{uj,0) ~ (— iu;)^^ proposed in 
Refs. 18 and 10 is valid only for C < 0. From Eq. (25) we obtain for the permittivity 



e(0,0) = l + 4Tre'^N{eF)d{0) 



- 1, 



0<C<2 
C>2 



(69) 



In the framework of the general analysis, e(0, 0) can diverge according to a law that is 
different from (obtained by cutting off the metaUic behavior of e(0, q) ~ on the scale 
q ~ More than that, e(0, 0) can even be finite in the limit ^ — > oo (see the discussion 
in Refs. 23, 40, and 41). 



5. BASIC STRUCTURE OF THE THEORY 

It is convenient to begin the construction of the theory by analyzing the localized phase, 
while the metaUic state will be obtained as a result of the instabihty of the locahzed phase. 
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5.1. Spectrum of the operator L in the localized phase 



Let M be the set of values of the index s that enumerates the eigenvalues of the 
operator L. We shall show that in the localized phase the decomposition 



M = Mo ® Ml e 



such that 



A., 



s e Mo 
s e Ml , 
s e Moo 



(70) 



(71) 



(Fig. 4) is valid. The set Mq is not empty, since it contains the element Aq ~ o;, related with 
the diffusion coefficient. We shall show that it is not the only element. According to Eq. 
(57), 0kk'(^) contains the singularity 1/uj, associated with two diffusion poles. In Eq. (38) 
this singularity originates from terms with s G Mq. Comparing these two representations 
and taking into account the fact that the diffusion pole at q — corresponds to the term 
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with s = in Eq. (38), we obtain 



l + d(k + k')(k + k')2 ^1^, l + ^,(q) ' ^ ^ 

where Mq is the set Mq without the clement s = 0. Since d{q) diverges as ^ oo [see Eq. 
(64)], the left-hand side of Eq. (72) contains a (5-function singularity at k + k' = 0, which 
terms /k/k' on the right-hand side of Eq. (72) with k' = — k cannot have at the arbitrary 
point k. The same is true for the sum of a finite number of such terms. The example of 
the Fourier expansion 

,i(k+k')-x 



l + d(k + k')(k + k')^ 



Y,A^e'^^+^>-'' (73) 



shows that the pole term in Eq. (72) can be reproduced by an infinite number of terms 
/k/k' and that this does not require a complete system of functions (eliminating from the 
sum in Eq. (73) terms with small x leads to the appearance of a smooth component, but 
it does not change the singularity at k-|- k' = 0). It is clear now, that the set Mq contains 
an infinite number of elements, but generally speaking it does not coincide with the set M. 

Sadovskii [10, 44] proposed a localization criterion according to which a nontrivial so- 
lution of the homogeneous Bethe-Solpeter equation appears in the hmit a; — > 0. We can 
make a stronger assertion: An infinite number of such solutions appears at the transition 
point. 

The following decomposition of the operator L follows from Eqs. (36) and (12): 

L — Lygg -\- Lsingi C"^^) 

L.epV'k = MV'k + AEk(q)Vk - ^ E ?7^k'(q)V^G'k(q)AG'k'(q)V'k', (75) 

k' 

Lsingi^v = E <'?'(q)VAG'k(q)AGk'(q)V'k'. (76) 
k' 

In the localized phase the diffusion pole in gives a 1/a; singularity, 

L = Lreg + —, (77) 

where the limit a; — >■ has been taken in the operator Li and terms of higher order in a; 
are included in Lreg- From Eqs. (76), (12), and (72) we obtain the following representation 



for Li 
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_ /■l"(q)v'^aW 

'^k+q/2'^k-q/2 

It is clear, that the eigenvectors of the operator Li, corresponding to nonzero eigenvalues, 
lie in the subspace constructed on the vectors 5'k^'*(q), and the number of eigenvalues is 
equal to the number of elements in Mq. The nonzero eigenvalues of Li correspond to the 
eigenvalues ~ l/cu of the operator Lgmg- 

The overall picture is as follows (Fig. 4). The operator Lging has an infinite number of 
eigenvalues ^ l/u and an infinite number of eigenvalues equal to zero. When the operator 
Lreg ~ 1 is added, the eigenvalues ~ l/tu change very little and form the set M^o of the 
operator L; the zero eigenvalues become, generally speaking, of order one, forming the set 
Ml, but infinite number of them remain ~ u and lie in the set Mq. The number of elements 
in Moo is equal to the number of elements in Mq, no assertions can be made with respect 
to the set Mi, but this is not important for what follows. 



5.2. Relation between L and Lsing 
We now introduce the spectral representation for the singular part of the operator 

Lsing — Li/u) 

f \ / I / seMo©Mi 

Li^2^ \us)Vs{us\, = I ^ 1 s e Moo ^ ^ 

and find a relation between L and Lf^ing, regarding L^eg as a perturbation. For s G Mqo 
the ordinary perturbation theory can be used, since all differences of the eigenvalues ~ l/cu 
and a regular expansion in powers of uj is obtained: 

I \ I \ , \^ \Us'\Lreq\Usl \ \ 

\es) = \us) >, \us'), 

P^, Vs - Vs' 

\s = — + {Us\Lreg\Us), S e Moo. (80) 

CO 

For s e Mo ® Ml we seek the eigenfunctions of L in the form 

|e)= E Cs\us) +u; J2 ^s\us), (81) 
seMoeMi seMoo 

where Cg, Dg ~ 1. Substituting (81) in the eigenvalue equation, we obtain a system of 
equations for Cg and D^, which can be solved by iterations in uj. Eliminating D^, we obtain 
to first order in uj 

{K^ss' - Tss')Cs' =0, s e Mo © Ml, 

s'eMoeMi 
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rp 1^, \j u, \ , , (-Ug I Lreg \Us"){Usii\ Lreg \Us') , . 

i.e., an ordinary secular equation taking into account the first correction from transitions 
into states witfi s e M^. 

5.3. Mechanism for satisfying the Ward identity 

We now demonstrate the cancellation of the singular contribution ~ l/o;, associated 
with the diffusion pole in C/kk'(*l)' the right-hand side of the Ward identity (13). The 
specific form of Lreg was not used in Sees. 5.1 and 5.2. To determine L^eg in the form (75) 
with M = we have 

1 Y: ^AG^{q)e^:\q) =0, se M^, (83) 

k 

since y^AG'k(q) is the exact eigenfunction of L belonging to the set Mq. By virtue of the 
relation (80), the difference of \us) from \es) for s G Moo is of order uj for any Lreg, whence 

1 E 7AGk(q)4^^(q) = 0(a;), « e Mo^. (84) 
-'^ k 

Comparing Eqs. (76) and (79), we have 

-[/-^^(q)VAGk(q)AG'kKq) = ;^ E 4^^(q)^.(q)4'^(q), (85) 

so that the singular contribution on the right-hand side of Eq. (13), taking into account 
Eq. (84), has the form 

T7L^kk'^G'k'(q) = -- E ^T^=T^^v^^k'(qK'(q) = — — (86) 

^ k' ^seM^ VAGk(q) k' ^ 

and the Xjuj singularity cancels. For the same reason, there will be no singularities on the 
right-hand side of Eq. (13) as the transition into the metallic phase is approached. In 
this case the spectrum of the operator L has the same structure (Fig. 4) with a; replaced 
by — the characteristic value of the diffusion coefficient (Sec. 6.2) — and we obtain 
0{Do)/Dq on the right-hand side of Eq. (86). 

5.4. Symmetry approach 
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The symmetry of the system is clearly expressed in the properties of the operator Lging: 

(a) The spatial uniformity in the mean makes it possible to introduce the three-momentum 
notations (Fig. Id) and to introduce the operator L in general and the operator Lgmg in 
particular. 

(b) The isotropy in the mean, combined with time-reversal invariance, guarantees that 
L and Lging are symmetrical and the existence of orthonormal bases of eigenvectors for 
them. 

(c) As a result of time-reversal invariance, Lging has a high symmetry, manifested in the 
existence of an infinite number of zero mode^. 

The decomposition (77) represents the operator L as a sum of the operator Lsing with 
a high degree of symmetry and a regular operator L^eg of a general form. It is similar 
to the decomposition (1) and is convenient for symmetry analysis. The condition on the 
transition point will be determined below and the origin of the parameter r will thereby 
be determined. 

Following Sec. 1, we consider the response of the system to a perturbation 6Lreg of a 
general form. Not all changes in the system will be important. We decompose the change 
in the operator L into two parts 

SL = 5Lx + 5L„ (87) 

where SL\ changes the eigenvalues and 6Le changes the eigenf unctions of L. For an infinites- 
imal 6L such a decomposition is trivial — 6Lx and 6Le are the diagonal and off-diagonal 
parts of the operator 6L in the representation of the eigenvectors le^). Changes of the SL^ 
type do not change the eigenvalues of L and therefore the diffusion coefficient D{uj,q), 
directly related with Ao(q). The diffusion coefficient determines uniquely the location of 
the system — in a localized phase, in a metallic state, or at the transition point. It is clear 
that the changes 6Le do not drive the system out of the transition point, they only displace 
the system along the critical surface [1]. Such displacements do not lead to nonanalyticity 
of the physical quantified and they can be ignored. The critical exponents obtained by 
motion along the normal to the critical surface are identical to the exponents obtained un- 
der an arbitrary nonzero angle to the tangent plane. Similarly, in perturbations of the 6Lx 
type, the part corresponding to a change in with s G Mi © M^o need not be considered. 

Only the changes in the eigenvalues from the set Mq are important, and their re- 
sponse to a perturbation is indeed nontrivial. Let the system lie deep in the localized phase. 
A small perturbation 6Lreg does not drive the system out of the state of localization and 
preserves the proportionahty As ~ a; for s G Mq. On the other hand, a perturbation SLreg 
of a general form possesses nonzero matrix elements with the respect to the eigenvectors |es) 
of the subspace Mq and should lead to small but nonvanishing, in the limit uj ^ 0, values of 



■^Their presence (see Sec. 5.1) is associated with the existence of a diffusion pole for k + k' ^ 0, which 
foUows from Eq. (9) (see Sec. 3). 

^Singularities associated with a change in the type of phase transition — for example, a second-order 
phase transition into a first-order phase transition — can occur on the critical surface. We assume that 
the system is far away from such singularities. 
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Xg. The resolution of this contradiction will lead to the self-consistency equation (Sec. 5.6). 



5.5. "Rotation" of the singular operator 

To formulate an adequate language for the further discussion, we shall consider the 
following problem of the "rotation" of a singular operator. 

Let the decomposition (77), where a; — > 0, be valid for the operator L. The operator 
Lreg acts in the space Jl, while the operator Li has nonzero eigenvalues ~ 1 in the subspace 
Qi, which is a part of Q, fl — Qq ® Qi. This justifies retaining in Eq. (77) two terms 
of different orders. Let SLi be a perturbation of the operator Li. If this perturbation 
is of a general form, then the addition SLi/uj to the operator L can be studied by the 
standard perturbation theory and gives corrections ~ l/cj. Let the perturbation SLi be 
such, however, that the operator Li + 5Li has the same properties as the initial operator 
Li. Then the dimension of the subspace Jli remains the same and only a "rotation" of the 
operator Li occurs; in this case SLi has no nonzero matrix elements in Qq- What is the 
result of such a perturbation in the subspace QqI 

Let fjs and \us) be the eigenvalues and eigenvectors of the initial operator Li. The 
operator Lreg can be neglected in the "upper" subspace fli, and in the "lower" subspace 
a secular equation in terms of the matrix elements {us\Lreg\us') should be written out. The 
perturbation 5Li produces the change 5us ~ 5Li of the eigenvectors \us) and the matrix 
of the secular equation is determined by the elements 

{Us\Lreg\Us') = {Ug + SUs\Lreg\Us' + SUg') = {Us\Lreg + SV\Us>). (88) 

The qualitative result is that a limitation of the form of the operator 5Li weakens its action 
on the lower subspace: The effective perturbation 6V appears to be ~ SLi instead of SLi/u 
for the operator of general form. 

The change in \us) in the subspace fii can be calculated by the standard perturbation 
theory, since all differences of the eigenvalues ~ 1 and a series in the small parameter arises: 

\us) = \us) + 2^ — z z \us'), sevti. (89) 



s 



% - Vs' 



An arbitrary choice of \us) that is compatible with the orthogonality relations can be made, 
in view of degeneracy, in the subspace Qq- To first order in 5Li we can set 

\us) = \us} - 2^ \us'}, s e Qo- (90) 



Vs 



Substituting the expression (90) into Eq. (88), we obtain for the matrix elements of the 
effective perturbation 



/- Ix-Cn- \ Sr^ {Us\SLi\Us"){Us"\Lreg\Us') + {Us'\SLi\Us"){Us"\Lreg\Us 

{Us\OV \Us') = — 2^ Z 

s"ef2i Vs" 



(91) 
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5.6. Self-consistency equation 



It is now easy to understand how to resolve the contradiction stated in Sec. 5.4. The 
perturbation SL^^g produces the change 5d{q) in the diffusion coefficient (62), which in view 
of the relation 

^ , iy(k,k',q)^k' 



2 



Nj^-iuj + D{oj, k + k') (k + k') 

d'^q W(k, -k + q, q)^-k+q _ -f/i'i/'k 



(92) 



{-iiu) J {2ttY 1 + d{q)q^ uj 

gives the following change in Li: 

Rotation of the subspace of the operator Lging produces in the subspace Mq the effective 
perturbation 5V , which in zeroth order in uj compensates SLreg- 

Introducing into Eq. (82) the small changes SLreg and SLi [the latter enters via the 
change in the eigenfunctions (90)], we obtain for the matrix of the secular equation 

T,s' - {Us\f + SLreg + SV\Us>) S, s' G Mq © Mi, (94) 

where the overbar denotes the unperturbed value, and T and SV are determined by the 
expressions (82) and (91) (with the substitution Qi M^o, ~^ © Mi). In the terms 
~ we confine ourselves to zeroth order in the increments. The choice of the vectors \us) 
in the subspace Mq ® Mi is arbitrary in view of the degeneracy. We choose them so as to 
diagonalize the matrix T — then, to zeroth order in cu, they are identical to the eigenvectors 
jcs) of the operator L [see Eq. (81)]. Since the eigenvalues of the matrix T are identical to 
the eigenvalues of L, we have 

T,s'^XAs' + {es\SLreg + SV\es'), s,s'eMo®Mi. (95) 

For infinitesimal 6Lreg and SV, the diagonal elements of the matrix T determine the eigen- 
values of the operator L 

Xs^Xs + {es\SLreg + SV\es), s e Mq ® Mi (96) 

and the off-diagonal elements determine the corrections to its eigenfunctions; the latter 
correspond to the perturbations of the type SLe (Sec. 5.4) and can be dropped. For 
constant je^) it is possible to switch in Eq. (96) from infinitesimals to finite increments. 
Further, the changes in Ag in the subspace Mi can be ignored (Sec. 5.4). Finally, we note 



25 



that fixing Ao(q) for all q means fixing the diffusion coefficient, which in turn determines 
all Xs{c[) with s e Mq, which can be reconstructed according to the binary decomposition 
(72). Therefore if Eq. (96) is satisfied for s = 



s"eMo, 



Vs" 



(Eqs. (55), (91), and (80) were employed), then it is automatically satisfied for all s from 
Mq. It is easy to show (see Appendix) that the expansion in q of the right-hand side of Eq. 
(97) contains only even powers of q, and terms ~ g° are absent in each of the two terms. 
Setting 

{eo\SLreg\eo) = -iq^5f{q) (98) 
and substituting the expression (93) into Eq. (97), we obtain 

q) - D(u, q) = 6f(q) + Q6d(q), (99) 

where Sf{q) and B{q,q) are regular functions of general form of the arguments q^ and g^. 
The quantity H^(k, k',q) in Eq. (92) can be expressed, by virtue of Eqs. (76), (12), and 
(59), in terms of functions which are regular at the transition point. This makes it possible 
to vary only d{q) on switching from Eq. (92) to Eq. (93). The equation (99) contains 
the diffusion coefficient on the right- and left-hand sides and replaces the self-consistency 
equation (16) of the VoUhardt-Wolfie theory. 



5.7. Condition on the transition point 

In deep of the localized phase D{uj,q) and D{uj,q) vanish at u; = and Eq. (99) 
determines the change 6d{q) = —Q~^Sf{q) for the prescribed perturbation 6Lreg- Making 
the small changes SLreg, we obtain the corresponding changes Sd{q), which preserve the 
proportionality D{u!, q) ~ {—icu). This situation remains as long as there exists an operator 
inverse to Q, i.e. as long as all eigenvalues of Q are nonzero. Let a nonzero eigenvalue of 
the operator Q appear at some point in the course of the motion from the interior of the 
localized phase. As we shall see below, such a point corresponds to the physical notions of 
the Anderson transition. 

The divergence of d{q) in the hmit ^ — > oo (Sec. 4) means [see Eq. (100)] that at the 
transition point the operator Q vanishes entirely or on some subspace. For this reason, it 
should be kept in mind in the analysis that many or even all eigenvalues of the operator 
Q can vanish simultaneously at the transition point. It is convenient to introduce the 
critical exponent Sg '>0 for each of them: 

QMq)^ f^sMl), f^s^r^' (101) 
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As d{q) — s> oo, the changes of the function B{q, q) cannot make the operator Q finite 
and therefore they do not drive the system out of the critical point; they only displace 
it along the critical surface and can be ignored. So the function B{q, q) is considered as 
independent of r. In this case, the equality 

/(0^^(^'^X^) = O (102) 

cannot be satisfied for any function (f){q). Indeed, it corresponds to the presence of a zero 
mode for the transposed operator (and therefore for the operator Q itself) not only at 
the transition point but also in a finite interval around it. 



6. SOLUTION OF THE SELF-CONSISTENCY EQUATION 

6.1. Claissification of the possible solutions 

The self-consistency equation for the metallic phase can be derived only by making 
specific assumptions about the functional form D{uj,q). For this reason, it is convenient to 
examine several cases which exhaust all possibilities. 

(a) Let there be among the exponents 5s in Eq. (101) a maximum exponent (for 
dcfinitcness, 5o)i i-e., among the set of soft modes, one mode is the softest. Then, as the 
transition is approached, the component const 0o(?); contained in Sf{q), will give rise to 
an anomalously large response constT~^°(j)o{q) in the function 6d{q). For this reason, near 
the transition the solution can be sought in the form 

D{uj, q) = DoIMq) + ^(q) « Mq)- (103) 

For 0o(?) to dominate for all values of q, it is necessary that 0o(O) 7^ 0, which we shall 
assume is the case. 

(b) Let several exponents have the maximum value 5o = Si = ... = 6p, and let at 
least one of the functions (j)o{q), 4>i{q), ■■■4>p{q) [for example 4>o{q)] different from zero 
at q — 0. Then near the transition 

D{uj, q) = Do[0o(g) + Ci0i(g) + ...+ 

+Cp</.p(g)+^(g)], (104) 

where Ci ~ C2 ~ ... ~ Cp ~ 1, (p{q) -C 4>o{(l)- 

(c) If for two eigenfunctions we have 0o(?) ~ (f^°: 0i(?) ~ 9^"^ in the limit g ^ 0, and 
no > ni. So > Si, then in the expansion of D(u, q) in (j)s{q) both functions must be retained. 
Although the coefficient of 0o(?) grows more rapidly near a transition, the function (f)i{q) 
dominates for small values of q. In the general case, d{q) must be sought in the form of the 
expansion (64) with arbitrary 
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Actually, as we shall see below, the case (b) is realized (Sec. 6.3). But the analysis 
of this case is virtually identical to the simpler case (a) (Sec. 6.2), which reproduces the 
solution of the self-consistent theory of localization [18]. Analysis of the case (c) requires a 
special mathematical apparatus (Sec. 7), and there is no need to extend this analysis for 
the metallic phase, since the solutions different from (b) do not exist. 



6.2. Case of a single dominant mode 

We seek the solution in the form (103). The definition of the operators Li and Q in 
Sec. 5 presumed that only the dependence on u in the localized phase is investigated. To 
investigate the dependence on uj and r it is necessary to take into account the fact that 
near the transition and in the metallic phase the magnitude of the diffusion denominator 
is determined by the parameter Dq ^ lu. Making the decomposition 



-iu + D{u;,q)q'^ Dq [(l)o{q)q'^ 0o(?) H^/-Do + 0o(?)?^ + </?(?)?' 
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(105) 



we write Lging in the form 



Lsing = n^^^ (^06) 

where Li and 6Li correspond to the first and second terms in the braces in Eq. (105). 
Substituting 6Li into Eq. (97) gives, instead of Eq. (99), the equation 

D{uJ, q) = Tf{q) + Qnifiq) - ^Quq'', (107) 

^ I rfl^r-^r •^/nl^r~^~2i ^(g)^ (108) 

J {2t:Y (t)Q{q)[-iuj / Dq + (t)Q{q)q^\ 

where we have neglected ip{q) in the denominator of Eq. (105), and have written 5f{q) in 
the form Tf{q). We also accept that the operator Li corresponds to the limit cu, r — > (see 
Fig. 5 below) and put D{uj,q) = 0. For the decomposition (103) to be unique, we require 
that ip{q) satisfy the condition 

expressing the requirement that (p{q) "not contain in itself the component const • (polq) 
{(poil) is an eigenfunction of that corresponds to the eigenvalue /Iq). Forming the scalar 
product of the expression (107) with 0o(?), we obtain 



Do{cl>0: 0o) = r(0o, /) - ^(00, QRq-'): (110) 
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where the last term is different from zero in view of the impossibihty of Eq. (102), and 
the term with ip{q) is absent since {4>o, Qr^p) = Qji4>o) = l^oi^^ 4>o) — 0- Written out in 
detail, Eq. (110) has the structure 

° Do J (27r)'^ MqWH^/Do + MqW] ' ^ ^ 

For d > 4 the integral is determined by large values of q, and —w/Dq in the denominator 
can be neglected. For d < 4 the integral is determined by small and we can set = in 
the slowly varying functions B{q) and 0o(?); the region of integration can be taken infinite 
and the integral can be made dimensionless. The result for both cases can be written in 
the unigue form 



Do^At + B[-—J , (112) 

introducing exponent v according to Eq. (18). The equation (112) has two types of 
solutions: in the metallic phase Dq — const 7^ as a; —> and Eq. (112) gives Dq — At 
in accordance with the value s = 1 for the conductivity exponent (18); in the dielectric 
phase Dq = {—iuj)C,'^ and ^ ~ t~'^ in accordance with the definition of the exponent of the 
localization length (in the case at hand, the configuration of exponents corresponds to 
the case Fig. 3a). The equation (112) and the values of the indices s and v are identical 
to those obtained in Ref. [18]. 

For the case d > i Eq. (112) reduces to a quadratic equation and it is easy to trace 
how the solutions are selected (Fig. 5). For u = the terms Dq = At and Dq = are 
intersected (Fig. 5a); for finite u), the degeneracy is removed by the amount ~ cj^/'^'^'^^ (Fig- 
5b), and of the two branches, only one satisfies the condition ReD{uj, q) > 0, following from 
Eqs. (31) and (35) and the non-negativity of {pePe+uj)<i- Choosing the indicated branch 
and passing to the limit c^; ^ 0, we obtain finiteness of Dq only on one side of the transition 
— for definiteness, for r > (Fig. 5c). 

From Eq. (107) we have for the function (p{q) 

V{q) = Q-r'Pi. (DoMq) - rfiq) + '^Qr^'") > (US) 

where P± is a projection operator onto the subspacc which is orthogonal to 0o(?)- Since 
Qr ~ 1 (sec below), we obtain ip(q) ~ max{|r|, c<j-'^/'-^'^+^^}, which justifies the assumption 
<f{q) ^ 0o(q')- For d > 2 the integral in Eq. (108) is determined by large vahics of q for 
any regular function 99(g) and all eigenvalues /if of the operator Qr are found to be of 
order unity. For the operator Q from Sec. 5 (which differs from Qr in the localized phase 
by the factor ^~^) this means that all Hs vanish according to the same law. Therefore the 
assumption that one mode predominates is not confirmed by the result and actually the 
case (b) of Sec. 6.1 is realized. 



6.3. Case of several dominant modes 
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Figure 5: Selection of the solutions of the self-consistency equation: a — Intersection of 
the terms for c<j = 0; b — splitting into physical and unphysical branches for a; > 0; c — 
behavior of Dq{t) for the physical branch at a; = 0. 
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We seek D{uj,q) in the form (104), where the choice of the function (p{q) is fixed by 
the conditions (0o;<^) = 0, Cj = const(r) as r — > (if </? is required to be orthogonal to 
4>i,---,4>p, then the coefficients Ci are functions of r, and this leads to inconveniences in 
defining the operator Li corresponding to the limit uj,t ^ and not depending on r). 
Using instead of 0o the "correct" linear combination 0o + Ci0i + ••• + Cp4'p and repeating 
the agrumcnts of Sec. 6.2., we arrive at equations of the type (107) and (108); forming the 
scalar product of the first equation with 4>q, we arrive at Eq. (112) with all consequences 
following from this. Once again, all eigenvalues of Q vanish according to the same law and 
the limit p — oo must be taken in Eq. (104), i.e., all (psio.) must be inchided in the correct 
linear combination. Forming the scalar product of the analog of Eq. (107) with 0i, ^2, 
we obtain a system of equations for Q: 

DoCSi, (f)i) = r{4>i, f) + nf{4>i, ip)- 

-^{^i,QRq-'), i = l,2,..., (114) 

where //f ~ 1. The function ip{q) is found to be ~ r, and in the limit a; — > 0, it has a 
discontinuity at r = 0, i.e.. 

Substituting the expression (115) into Eq. (114), we obtain in the limit a; — >^ the equations 

= r(0„/) + r/iffif (0„(/),) +r'(0^,QR9"') (116) 

for the metallic and dielectric phases, respectively. For any Ci the equations (116) can be 
satisfied by appropriate choice of Bf"^ and B^^ , i.e., the coefficients of the correct linear 
combination are completely arbitrary. The meaning of this arbitrariness will be explained 
in Sec. 8. 

Finally, D{u!,q) near the transition has the form 

D{u;,q) = Dodiq), 

( r, r > cui/(2-+i) 

Dor^l cuV(2-+i), |r| ~ a;i/(2.+i) ^ (^60) 
[ (-io^)|r|-2^ -r>o^V(2-+i) 

where the function d{q) = d{q)/d{0) varies on the scale g ~ A. This result, obtained for 
D{(jj, q) defined as in Eq. (55), is also valid for the observed diffusion coefficient Dobsi^-, q)-, 
since the renormalizations associated with the functions -B(q) and C(q) in Eqs. (53) and 
(54) either contain no anomalous dispersion or they are small. 
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7. UNIQUENESS OF THE SOLUTION 

In this section the self-consistency equation (99) in the locahzed phase is investigated 
assuming for d{q) an expansion of the general form (64). 



7.1. Method of supporting points 

In what follows, integrals of the form 

{2-nY + ^"'g^ + ... + _^ ' 



an > (117) 



are essenrial. The asymptotic expressions of these integrals in the limit ^ oo are cal- 
culated by the method of " supporting points" . We choose an appropriate scaling in order 
to make the integrals dimensionless, making the substitution q = $,~H, and removing the 
common factor from the denominator; as a result, the indices a<j become as — a — 2sb. 
By choosing appropriate values of a and b the exponents in two terms of the denominator 
in Eq. (117) can be made to be zero and the remaining exponents become negative. Then 

/Jdj. j.2k 
Jll ^ ^-a-b(d+2k) . ON 

(27r)'^t2si +^2S2 + y: ^^s-a~2sbf2s ■ K J 

S^Sl,S2 

The integer numbers si and S2 should satisfy the condition 2si < d+2k < 2s2, guaranteeing 
that the integral converges after the sum over s is dropped; to avoid uncertainties, we 
assume that d is noninteger, and pass to the limit of integer d in the final results. The 
procedure described above admits a simple geometric interpretatioijf]. Let us plot a sequence 
a„ against n (Fig. 6a), mark on the abscissa axis the point = {d + 2fc)/2, and construct 
"an upper tangent" at the point Xk to the set of points (n, a„). If we imagine, that the 
points are represented by nails, then this construction is made with the aid of "stick" (solid 
line in Fig. 6a) and "rope" (dashed line). The numbers of the points, on which the upper 
tangent "lies", determine Si and S2, and its equation a = a + 2bn determines the parameters 
a and b. 

Constructing the broken line, consisting of segments of the upper tangents (Fig. 6b), 
we obtain a "convex envelope" a{x), in terms of which the result (118) assumes the form 

4~r"^-\ xk = ^. (119) 

By construction the function a{x) is increasing and convex (in the not strict sense). For a 
bounded sequence with maximum at n = no, a{x) is strictly increasing for x < hq and 



'Similar constructions arise in the investigation of Burgers equation [45]. 
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Figure 6: a — Construction of the upper tangent at the point Xk to the set of points 
{n, an); b — convex envelope a{x) for the sequence «„. 

constant for x > no (for Xk > no the supporting point S2 hes at infinity). For a strictly 
increasing and strictly convex sequence the following inequalities follow from Eq. (119): 



/o > i^i > ^2 > ... » 4 » 4+1 > 
4i42 ^ 4i-i42+i hi-2lk2+2 <c ... , ki < k2. 



(120) 



'i2r 



For an arbitrary sequence some of the strong inequalities are replaced by weak inequali- 
ties. In what follows, for definiteness, we proceed from the strong inequalities in Eqs. (120) 
and (121), having in mind that the results remain valid in the order of magnitude for the 
general case. 



7.2. Symmetrization of the operator Q 

We set in Eq. (99) Sf{q) = STf{q) and expanding all functions in series 



2k ~2k' 



f{q) = E hQ'", Hq) = E ^dkq'", B{q, q) = E Bkk'q^'q 



k=0 



we obtain in the limit a; — > 



fe=0 



k,k'=0 



(122) 



where 



-^^fk — E ^l^k' E ^k'+k"Sdk", 
k'=0 k"=0 



I 



d'^q 



-,2fe+2 



(27r)'^[l + d(?)?2]2- 
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(123) 



(124) 



The matrix B = \ \Bkk> \ \ and the column matrix / = | | are of general form with elements 

~ 1. The column matrix / = B^^f has the same properties. Multiplying Eq. (123) by 
B'^ we obtain 



-St 



fo 
fi 

h 



h 
h 



h 

h 
h 



h 

h 
h 



Sdi 
Sd2 



(125) 



i.e., an equation of the type (99) but with the symmetrized matrix of the operator Q. 

Using the expansion (64) for d{q), the integrals (124) acquire the form (117) with the 
exponents 

an = max{k,}i(3ki + fik2)k^+k2=n (126) 

and an extra in the numerator. The exponents j3n are nonnegative and increase not 
faster than 2n (Sec. 4). This guarantees the condition q;„ > and make is possible to 
construct a convex envelope. 



7.3. Inversion of the operator Q 

Restricting the upper limit of the summation in Eq. (122) by some finite n, we obtain 
in Eq. (125) a system of equations of finite order that can be solved by Cramer's rule. The 
determinant of the matrix Q in Eq. (125) consists of all possible products of the form 

hohi+ih2+2---h„+n, (127) 

where fco, fci, A;„ is a permutation of 0,1,. ..,n. We separate in Eq. (125) the pair 
Iks+sh'^+s' with s < s'. If ks > kgi-i then it follows from Eqs. (120) and (121) that 

h'^+shs+s' ^ hs+sh^i+s' (128) 

and the product (127) can be increased by interchanging kg and kg', without touching the 
other fcj. It is obvious that in the maximum product among the products (127), which 
determines the order of magnitude of the determinant Q, should have ko < ki < ... < kn, 
whence ko = 0, ki = 1, kn = n and therefore 

detg ~ lohh-hn- (129) 

The minor of the matrix Q, obtained by crossing out the ith row and the jth column, 
consists of all possible products of the form 

4o4i+i---4j_i+(i-i)4,+i+(j+i)---4„+n, (130) 

where /cq, A-'i, kj^i, fcj+i, A;„ is a permutation of 0, 1, i— 1, i + 1, n. In the maximum 
product these two sequences are identical. It is easily verified that 

(131) 
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Solving Eq. (125) by Cramer's rule and using Eqs. (129) and (131), we obtain 

ddk-^dT- — , A; = 0, l,...,n (132) 

and, using Eq. (119), the result can be expressed in terms of the sequence a^. For a convex 
sequence f3k-, we have from Eq. (126) = 2/9m,, «2m+i = Pm + f^m+i, which can be written 
in the form ak = 2/9^/2, if the sequence (3k is additionally defined at half-integral points by 
the relation /3k+i/2 = {Pk + Pk+i)/2. Since the values of /3k+i/2 he on the convex envelope 
(3{x), for arbitrary x we obtain 

a{x) = 2P{x/2). (133) 

This result remains valid for an arbitrary sequence Pk- To prove this, it is necessary to 
introduce an auxiliary convex sequence Pk = P{k) > Pk and note that replacing Pk by Pk 
does not change the value of the integrals Ik- For convex Pk we have 

Pik + ip)^il-ip)Pk + ipPk+i, 0<ip<l, (134) 

which makes it possible to switch from the convex envelope directly to the values of Pk- In 
the general case Eq. (134) is correct with Pk replaced by Pk- Setting 

d = Am + 4-0, m — integer, < -0 < 1 (135) 

we obtain from Eqs. (132), (119), (133) and (134) in the limit n — > oo 

5dk - 5re^(*^), 

gri^\ = / ~ 2?/;)/3^+fc + 2i)p^+k+i + Poo, < < | .^^gx 

\ Pm+k+l + Poo, 2 — — 1 

where the limit Poo — linife_*oo Pk is assumed to be finite in accordance with the considera- 
tion of the next section. 



7.4. Impossibility of unbounded growth of /3k 

For an unbounded sequence Pk the convex envelope Pk is strictly increasing and the 
hierarchy (120) continues to infinity. By virtue of Eq. (132), this means that Sdk diverges 
as n — > oo. To clarify the reasons for the divergence, we note that the off-diagonal part 
of the matrix Q in Eq. (125) under the conditions (120) and (121) can be regarded as a 
perturbation. Its eigenvalues in leading order are equal to l2k and they bunch up near zero 
in the limit /c — > oo. In the proof of Fredholm's theorem it is shown [46] that when the 
expansions (122) are truncated at the nth term, the (n + 1) maximum eigenvalues of the 
operator Q are reproduced; in the limit n — > oo, arbitrary small eigenvalues are reproduced 
and the response of the system to a small perturbation diverges. This situation occurs not 
only at the transition point but also in a neighborhood of the transition point (as long as 
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^ 3> A^^); it is unphysical, since the system is unstable with respect to an infinitely small 
perturbation of a general form. 

This result has important qualitative consequences, since it excludes the cases corre- 
sponding to Fig. 3b and proves the vahdity of the results D{u;, 0) ~ (— io;)^^ and e(0, 0) ~ 
in the localized phase. 



7.5. Change in d(q) as the transition is approached 

Expanding the numerator in Eq. (92) in powers of g^, we obtain integrals that can be 
calculated by the method of supporting points and which are of order ^-Pi^o)^^-P{xi)^ 
so on. We set 

where the terms C^^^k arise from the higher order terms in the expansion in (p and from 
corrections to the main scaling in the method of supporting points. We have changed the 
definition of Li in comparison with (77), in order to make it possible to separate the main 

singularity as r — > and introduce the operator Li, corresponding to the limit a;,T — > 0. 

In Eq. (137) it was assumed that (3^ are constant. Now, let the change 5t in the 
parameter r generate the changes 5^ and 5 (5k in the quantities ^ and (5k. Then 



Li + {i+ d[,f''-"\-i)dL, + (( + dO-"'h 



(138) 



where only the term with the minimum index yi is retained, and 5Li is determined by the 
expression (93) with 5d{q) of the form 

<^c?(g) = i:rt'''=+^lnm+i. (139) 

k=0 

Using as SLi in Eq. (97) the quantity Li(^ + 5^) — Li{^) we obtain instead of Eq. (99) 

5D{u;, q) = 5Tf{q) + C'^-'S^Riq) + ^^'''^QSdiq), (140) 

where 6D{uj, g) ^ as c<j — ^ 0. In the case of exact scaling, when 6d{q) = 0, the first two 
terms on the right-hand side cancel one another, whence yi = l/u. In the general case, 
they are of the same order of magnitude. Inverting the operator Q according to Sec. 7.3 
and comparing with Eq. (139), we have 

Sdk ~ St^-^^'='>^+^^''^ ~ e^'=+ilne%+i, (141) 
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and, consequently 



(142) 



where it is convenient to write 7^ in the form 



Ik = < 



< ^ < 1) 



iPk - h) + (^00 - /3(xo)) + (1 - 2^)(^„+fc_i - h)+ (0 < ^ < 



2' 



(A - + (^00 - ^k) + (1 - 2i,){^k-i 



m > 1) 

(0 < ^ < I, 
m = 0) 



(143) 



All combinations in parentheses are non-negative and 7^ > 0. For fixed ^ the definition 
of the exponents (3k in Eq. (64) is not unique: the coefficient of q^^ can be written as 
Ck^^^ and a change in (ik is equivalent to a change in Ck- The specific configuration of 
the exponents (3k make sense only if it remains unchanged when r changes. According to 
Eq. (142), for 7^ > a large change in the exponents occurs. Only if 7^ = 0, the changes 
~ 5t can be included in the changes of Ck- The condition 7^ = requires that 
all combinations in parentheses of Eq. (142) vanish and it fixes the only configuration of 
exponents, which is different for d > 2 and d < 2: 



Pi =(32 = (33 

(3o = (3i = (32 



d>2. 



d<2. 



(144) 

By definition, (3q = and for ci < 2 all exponents are equal to zero. This means that d{q) 
does not diverge and the localized phase remains for all r [21]. For d > 2, all indices can 
be made equal to 2, in accordance with the requirement Pi = 2 (Fig. 3a), by redefining ^. 
All eigenvalues of Q vary according to the same law and we return to be case (b) of Sec. 6.1. 



8. CHANGE IN SYMMETRY AT THE ANDERSON TRANSITION 

A change in Lreg gives rise to a rotation of the subspace Mqo of the operator I/^mg- This 
is analogous to a rotation of the magnetization vector M accompanying a change in the 
magnetic field H in a ferro magnet. This analogy is formalized in the form of Table 1. We 
shall give some explanation^. 

The operators Lreg and Lging have many degrees of freedom, many of which do not 
appear in the self-consistency equation. The important degrees of freedom are determined 
by the functions /(g) and d{q) [see Eq. (116a)], whose expansion coefficients 

f{q) = l + fiq' + ... + fnq''' + ..., 



similar, but not identical, analogy was discussed in Ref. [41]. 
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d{q) = l + d,q^ + ... + dnq^'' + ... (145) 

can be regarded as components of the unit vectors H and M. In the localized phase, small 
changes in them are related with each other by the operator Q, whose inverse is analogous 
to the magnetic susceptibility tensor Xij- 

The finiteness of the frequency a; smears the transition, similarly to the finiteness of 
the magnetic field in a ferromagnet. In the localized phase Dq u and in the metallic 
phase Do = const(cj), which is analogous to the appearance of spontaneous magnetization, 
i.e., the quantities Dq and u are analogous to |M| and |H|. In view of the qualitative 
character of the analogy, this identification is not unique. For example, any monotonic 
function F(|M|), equal to zero for |M| = 0, can be taken as the analog of Do; as a function 
of the magnetization itself, it has the form F{M'^), since a scalar must be formed from a 
vector. Finally, for small M^, it can be expanded in a series, and we obtain an analogy of 
Dq to M^. Similarly, is the natural analog for lj. 



TABLE 1. Analogy between a ferromagnet and a disordered system. 



Ferromagnet 


Disordered system 


Orientation of the magnetic field H: 

Components of the unit vector 
Orientation of the magnetization M: 
Components of the unit vector 
Squared modulus of the field 
Squared modulus of the 
magnetization 
Magnetic susceptibihty tensor Xij 
Paramagnetic phase 
Ferromagnetic phase 
Curie point 
T-T, 


Operator Lreg- 
Coefficients 
Space Moo of the operator Lsing 
Coefficients d„ 
Frequency uj 
Diffusion coefficient Dq 

Operator 
Localized phase 

Metallic phase 
Point of the Anderson transition 
Distance to the transition r 



In the analogy found, it is important that (a) the number of components of the vector 
M is infinite, since the number of expansion coefficients dn is infinite, and (b) the ferro- 
magnet is isotropic. The latter property is obvious from the fact that all eigenvalues of the 
"susceptibility tensor" Q^^ diverge at the transition point according to the same law and 
for small changes in dn and /„ they can be made equal by a r-independent linear trans- 
formation. The analog of collinearity of M and H in an isotropic ferromagnet exists with 
the following stipulation. In the case of the Anderson transition the "vector H" and the 
"vector M" fie in different subspaces and there is no natural method for estabfishing the 
mutual orientation of these subspaces. For this reason, for a fixed function f{q) the choice 
of d{q) is arbitrary (Sec. 6.3), in accordance with the arbitrariness in the choice of bases 
in the two subspaces. For the special choice Bf^ = in Eq. (115) we have /(g) = d{q), 
which corresponds to the choice of the "correct" mutual orientation of the bases. 
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The model of an isotropic ferromagnet with the number of components n — oo is well 
known in the theory of phase transitions and is the basis for the 1/n expansion [1]. Its 
critical exponents are known exactly. Specifically, for the exponents of the magnetization 
M and of the correlation length ^ ~ t~'^ we have 



/3 = 1/2 , d>2- 



u = < 



d-2 



1 
2 



2 < d < 4 



> 4 



(146) 



which, since s = 2/3, corresponds exactly to Eq. (18) 



9. CONCLUSIONS 

The approach to the theory of localization based on the formalism of a-models [22, 47, 
48] is currently considered to be the most rigorous approach. However, its rigor should not 
be overestimated. First, the degree to which the approximations employed in the derivation 
of the cr-models retain exact time-reversal invariance of the initial disordered system and the 
satisfaction of the Ward identity (13), was never studied. But these properties are vitally 
important for reproducing the pole structure of t/kk'(^)- Second, to take into account 
the spatial dispersion of D{uj,q), it is necessary to introduce into the Lagrangian of the 
cr-model additional gradient vertices, which grow anomalously at the initial stage of the 
renormalization group transformations [49]. The analog of such growth can be obtained 
from Eqs. (142) and (143), assuming for the initial configuration of exponents jSi = 2, f3k = 
{k> 2): 

Ae^-^N A; = 1,2,.... (147) 



dr In^ 

Growth of Pk with k > 2 indicates intensification of spatial dispersion of D{uj,q) as the 
transition is approached. In the language of the magnetic analogy (Sec. 8), it corresponds 
to a transformation of a uniaxial ferromagnet into an isotropic ferromagnet0. One can 
speculate, that the renormalization group transformations can transfer analogously the 
zero-component a-model into an infinite component model. These difficulties apparently 
are not important in low orders in e = d — 2, since for small e the Anderson transition falls in 
the region of weak disorder, for which the derivation of the cr-model is indeed substantiated. 

We now discuss the possible reasons for the disagreement of Eq. (18) with the result of 
Ref. 28 for the exponent s. The result of Ref. 28 for the permittivity e(0, 0) ~ ^ corresponds 
to the case C = 1 of the Sec. 4. This is also indicated by the expression given in Ref. 28 for 
the function A{r) from Eq. (66). For ( = 1 the exponents jSk increase linearly with k, and 
in accordance with Sec. 7.4, the instability with respect to an infinitesimal perturbation 



detailed investigation of the evolution of f3k requires a knowledge of the proportionality coefficients 
in Eq. (64) and (142). 
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of a general form arises. If the results of Ref. 28 correspond to the exact solution of some 
idealized model, then this model is unphysical. More likely, the approximations employed 
in the derivation of the a-model and the selection of diagrams in Ref. 28 destroy the pole 
structure of f/kk'(<l)- The fact, that results are the same for models with and without the 
time-reversal invariance, looks suspicious in this sense. Finally, in the derivation of the a- 
model for a large d it is necessary to introduce an artificial construction of weakly coupled 
granules. For this construction, because of the presence of artificial small parameters, the 
critical region can narrow anomalously and, as a result of the approximations, contract 
into a point. The results of Ref. 28 could correspond to some intermediate asymptotic 
behavior. 

There are some objections to the arguments of Efetov [43] that the diagrammatic ap- 
proach in principle does not "feel" the noncompactness, which in his opinion determines 
the main difference between the theory of disordered systems and the theory of phase tran- 
sitions. One can agree with the last assertion: Noncompactness is a consequence of the 
infinitesmall additions ±i5, determining the type of Green's function, which lead to non- 
perturbative contributions giving rise to the difference between the two indicated theories 
[17]. However, the nonperturbative contributions can be obtained from the diagrammatic 
technique [16]. The additions ±i5 play an important role in the separation of the diffusion 
poles, since as a result of these additions, the integration contour in Eq. (58) is confined 
between the poles of two Green's functions. 
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APPENDIX 

Expansion of the self-consistency equation in powers of q 

If we have for the operator L 

L(q)V^k = -^X^^kk'(q)V^k', 
k' 

^kk'(q) = ^k'k(q) = i>-k,-k'(-q) (^i) 

then it is easy to prove that: (a) the eigenvalues of L are even as a function of q, As(q) = 
As(— q); (b) the eigenfunctions e^\q) can be chosen so that e^'*^(q) = e^^(— q); and, (c) if 
there are several operators of the type (Al), then the matrix elements of one operator with 
respect to the eigenfunctions of the other operator are even functions of q. The operators 
L, Lreg, and Li are of the form (Al) and, by virtue of (-c), the right-hand side of Eq. (97) 
is even in q and can be expanded in powers of q^. 
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Since the operator Li in Sec. 5.1 is independent of lv, we obtain that in Eq. (86) 
0{u;) = and the definitions (75) and (76) are equivalent to the following definitions: 

Lregipi. = (ek+q/2 - ek-q/2)V'k + ^ E ^^kk' ( q) [^^^k' ( q) V'k - ^AG'k(q) AGk' (q) V'k'] , {A2) 

k' 

4.n,^k = ^E^kL''(q)[^Gk'(q)^k- V^Gk(q)AGk'(q)^k']. (A3) 
-'^ k' 

For the definition (A2) we have, on account of Eq. (39), 

Lreg\eo) = L^eg{const^AG'k(q) + 0{q)} = 

= const (ek+q/2 - ek--q/2)\/AG'k(q) + 0(q) = 0(q) 
and the contribution 0{q^) is absent in each term on the right-hand side of Eq. (97). 

References 

[1] Shang-Keng Ma, Modern Theory of Critical Phenomena, W. A. Benjamin, Reading, 
MA, 1976. 

[2] L. D. Landau and E. M. Lifshitz, Statistical Physics, Pergamon Press, N.Y. [Russian 
edition, Nauka, Moscow, 1976]. 

[3] A. Z. Patashinskii and V. L. Pokrovskii, Fluctuation Theory of Phase Transitions, 
Pergamon Press, N. Y., 1979. 

[4] A. A. Belavin, A. M. Polyakov, and A. B. Zamolodchikov, Nucl. Phys. B 241, 333 
(1984). 

[5] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 
[6] D. J. Thouless, Phys. Rep. 13, 92 (1974). 

[7] N. F. Mott and E. A. Davis, Electronic Processes in Non-Crystalline Materials, Claren- 
don Press, Oxford, 1971. 

[8] A. L. Efros, Usp. Fiz. Nauk 126, 41 (1978) [Sov. Phys. Usp. 21, 746 (1978)]. 

[9] M. V. Sadovskii, Usp. Fiz. Nauk 133, 223 (1981) [Sov. Phys. Usp. 24, 96 (1981)]; D. 
Belitz and S. R. Kirkpatrick, Rev. Mod. Phys. 66, 261 (1994). 

[10] M. V. Sadovskn, Sov. Sci. Rev. A. Phys. 7, 1 (1986). 



41 



[11] A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinski, Methods of Quantum Field 
Theory in Statistical Physics, Prentice-Hall, Inc., Englewood Cliffs, N. J., 1963. 

[12] S. F. Edwards, Phil. Mag. 3, 1020 (1958). 

[13] 1. M. Suslov, Zh. Eksp. Teor. Fiz. 83, 1079 (1982) [Sov. Phys. JETP 56, 612 (1982)]; 
84, 1792 (1983) [Sov. Phys. JETP 57, 1044 (1983)]. 

[14] I. M. Suslov, Zh. Eksp. Teor. Fiz. 92, 1433 (1987) [Sov. Phys. JETP 65, 806 (1987)]. 

[15] A. Aharony and Y. Imry, J. Phys. C 10, L487 (1977). 

[16] 1. M. Suslov, Zh. Eksp. Teor. Fiz. 102, 1951 (1992) [Sov. Phys. JETP 752, 1049 
(1992)]. 

[17] I. M. Suslov, Zh. Eksp. Teor. Fiz. 106, 560 (1994) [JETP 79, 307 (1994)]. 

[18] D. VoUhardt and P. Wolfle, Phys. Rev. B 22, 4666 (1980); Phys. Rev. Lett. 48, 
699 (1982); in Modern Problems in Condensed Matter Sciences, edited by V. M. 
Agranovich and A. A. Maradudin, North-Holland, Amsterdam, 1992, Vol. 32. 

[19] A. V. Myasnikov and M. V. Sadovskii, Fiz. Tverd. Tela (Leningrad) 24, 3569 (1982) 
[Sov. Phys. Solid State 24, 2033 (1982)]. 

[20] L. P. Gor'kov, A. 1. Larkin, and D. E. Khmel'nitskii, JETP Lett. 30, 228 (1979). 

[21] E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishman, Phys. 
Rev. Lett. 42, 673 (1979). 

[22] K. B. Efetov, A. 1. Larkin, and D. E. Khmel'nitskii, Zh. Eksp. Teor. Fiz. 79, 1120 
(1980) [Sov. Phys. JETP 52, 568 (1980)]. 

[23] K. B. Efetov, Zh. Eksp. Teor. Fiz. 88, 1032 (1985) [Sov. Phys. JETP 61, 606 (1985)]. 

[24] V. L. Berezinskii and L. P. Gor'kov, Zh. Eksp. Teor. Fiz. 77, 2498 (1979) [Sov. Phys. 
JETP 50, 1209 (1979)]. 

[25] H. Kunz and R. Souillard, J. de Phys. Lett. 44, L506 (1983). 

[26] S. Hikami, Phys. Rev. B 24, 2671 (1981). 

[27] H. Kunz and R. Souillard, J. de Phys. Lett. 44, L411 (1983). 

[28] K. B. Efetov, Zh. Eksp. Teor. Fiz. 94, 357 (1988) [Sov. Phys. JETP 67, 199 (1988)]. 
[29] B. Shapiro, Phys. Rev. Lett. 50, 747 (1983). 

[30] A. MacKinnon and B. Kramer, Phys. Rev. Lett. 47, 1546 (1981). 



42 



[31] M. Shreiber and B. Kramer, in Anderson Localization, Springer Proc. in Physics 28, 
edited by T. Ando and H. Pukuyama, Springer, Berlin, 1988, p. 92. 

[32] P. Wegner, Nucl. Phys. B 316, 663 (1989). 

[33] Y. Imry, Y. Gefen, and D. Bergman, Phys. Rev. B 26, 3436 (1982). 

[34] E. Abrahams and P. A. Lee, Phys. Rev. B 33, 683 (1986). 

[35] J. T. Chalker and I. H. Nahm, J. Phys. Condens. Matter 1, 3615 (1989). 

[36] L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media, Pergamon 
Press, N. Y. [Russian edition, Nauka, Moscow, 1982]. 

[37] E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, Part 2, Pergamon Press, N. 
Y., 1980. 

[38] D. Pines and P. Nozieres, Theory of Quantum Liquids, W. A. Benjamin, N. Y., 1966. 

[39] J. Ziman, Elements of Advanced Quantum Theory, Cambridge University Press, N. 
Y., 1975 [Russian translation, Mir, Moscow, 1971, p. 64]. 

[40] E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics, Pergamon Press, N. Y., 1981 
[Russian edition, Nauka, Moscow, 1979, Sec. 78]. 

[41] B. L. Al'tshuler, A. G. Aronov, D. E. Khmcl'nitskii, and A. I. Larkin in Quantum 
Theory of Sohds, edited by I. M. Lifshitz, Mir Pubhshers, Moscow, 1982. 

[42] I. M. Suslov, JETP Lett. 43, 704 (1986). 

[43] K. B. Efetov, Zh. Eksp. Teor. Fiz. 92, 638 (1987) [Sov. Phys. JETP 65, 360 (1987)]. 

[44] M. V. Sadovskii, Zh. Eksp. Teor. Fiz. 83, 1418 (1982) [Sov. Phys. JETP 56, 816 
(1982)]. 

[45] Ya. G. Sinai, Commun. Math. Phys. 148, 601 (1992). 

[46] V. S. Vladimirov, Equations of Mathematical Physics, M. Dekker, N. Y., 1971 [Russian 
edition, Nauka, Moscow, 1967, pp. 238-244]. 

[47] L. Shafer and F. Wegner, Z. Phys. B 38, 113 (1980). 

[48] K. B. Efetov, Adv. Phys. 32, 53 (1983). 

[49] V. E. Kravtsov, I. V. Lerner, and V. I. Yudson, Zh. Eksp. Teor. Fiz. 94, 255 (1988) 
[Sov. Phys. JETP 67, 1441 (1988)]; F. Wegner, Z. Phys. B 78, 33 (1990). 



43 



